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We review the Kauzmann paradox and what 
it implies about the configuration space en- 
ergy hypersurface for "structural glassform- 
ers." With this background, we then show 
how the relaxation expression of Adam and 
Gibbs qualitatively accounts for most of the 
phenomenology of liquid and polymeric 
glassformers including the strong/fragile 
liquid pattern, and the behavior of non- 
ergodic systems. Extended temperature 
range relaxation studies are consistent with 
a relaxation time pre-exponent on the quasi- 
lattice vibration time scale. When this 
boundary condition is imposed on Vogel- 
Fulcher-Tammann fittings, correspondence 
of T Q with T K is found for liquids with T g 
ranging over 1000 K. When it is imposed 



on the WLF equation C\ is obliged to be- 
come ~16, and the corresponding force-fit- 
ted C 2 provides a measure of the polymer 
fragility which is generally not available 
from thermodynamic studies. Systems 
which exhibit discontinuous changes in 
configuration al entropy on temperature in- 
crease, which include unfolding proteins, 
are briefly reviewed. 
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1. Introduction 

Walter Kauzmann wrote only one paper on the sub- 
ject of supercooling liquids and glasses, but it has 
proved to be one of the most influential papers in all of 
glass science [1]. There are many ways of presenting 
the main point he sought to make, but none are more 
telling than the graphical presentation chosen by Kauz- 
mann himself. Thus we reproduce in Fig. 1 the key 
figure from this seminal paper. We find this presenta- 
tion particularly attractive because it combines the key 
point which Kauzmann wanted to make, with a demon- 
stration of the concept of fragility [2] to which we have 
given much attention in our own work of the past 
decade. We will briefly review the paradox presented 
by the data contained in Fig. 1, and then devote the 
remainder of this article to the manner in which the 
approach of Gibbs and his coworkers to this intriguing 
problem leads to an effective (though so far inexact) 
resolution of the paradox, and at the same time, to a 
very broad qualitative understanding of the phe- 
nomenology of viscous liquids and the glasses which 
form from them. 



2. The Kauzmann Paradox and the 
Potential Energy Hypersurface for 
Glassformers 

In Fig. 1, the difference in entropy between the crys- 
tal and liquid at the melting point, A5V, is used as a 
scaling parameter to permit the simultaneous compara- 
tive display of the manner in which the difference in 
entropy between the liquid and crystal states, for six 
different substances, varies during their supercooling. 
All cases show positive slopes, reflecting simply the 
fact that liquids have higher heat capacities than the 
corresponding crystals. What is interesting is the mag- 
nitude of the slope. 

In the case of boron trioxide, now much spoken of as 
a "strong" liquid [2a,2c,3,4], the slope is such that the 
excess entropy of the liquid over crystal is only tending 
to disappear in the vicinity of K — which raises no 
concern at all. On the other hand, to various degrees, the 
other liquids in the figure show provocative behavior. 
In the case of lactic acid, which we would now call the 
most "fragile" of the six, the excess heat capacity of the 
liquid over crystal is causing the excess entropy to 
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Fig. 1. Kauzmann's presentation of the entropy crisis which bears 
his name. The figure shows the rate at which the difference in entropy 
between liquid and crystal, normalized at the fusion point, disappears 
at T is lowered towards absolute zero. For B 2 3 , now known as a 
"strong liquid," the liquid would always be of higher entropy than the 
crystal, even if the glass transition did not intervene at high T/T m , to 
change the heat capacity. At the other extreme, lactic acid loses its 
excess entropy so rapidly on cooling that if T $ did not intervene to 
arrest the loss, liquid would arrive at the same entropy as the crystal 
at 2/3 of the melting point. (This is the temperature usually associated 
with the temperature of the glass transition itself (the 2/3 rule which 
this set of data only weakly support).) Lactic acid is an example of a 
"fragile" liquid. Other examples of these plots for fragile liquids are 
given in Ref. [2]. 



decrease so quickly that a simple extension of the be- 
havior would lead that excess to vanish at a temperature 
which is only —2/3 of the fusion temperature. As far as 
can be told from the data, all that prevents this at first 
sight mind-boggling thermodynamic inversion from 
occurring, is the occurrence of a kinetic phenomenon, 
the glass transition at the temperature T s . At T s , the 
increasingly sluggish motion of the particles prevents 
further configurational changes to more ordered states 
from occurring as cooling proceeds at a fixed rate, and 
thus prevents equilibrium from being attained at lower 
temperatures. The consequent decrease in heat capac- 
ity, which is the signature of the glass transition, means 
the rate of entropy loss from the liquid is decreased, 
and the crisis is averted, as indicated in the diagram. 
The avoiding of a thermodynamic crisis by intercession 
of a purely kinetic phenomenon constitutes a paradox, 
the Kauzmann paradox, which must be resolved before 
glass science can be considered to have reached the 
state of a "mature field" of research. 



Two quite profound theoretical problems are pre- 
sented by the data of Fig. 1 . The first is the problem of 
constructing an equilibrium theory for the liquid state 
that contains an explanation of how, on infinite time 
scales, the system evolves so as to undergo a rather 
abrupt, if not singular, change in heat capacity at some 
temperature between the glass transition temperature 
and absolute zero. Part of this problem involves the 
interpretation of fragility of liquids and the coupling of 
vibrational to configurational degrees of freedom. The 
second is the problem of constructing a theory which 
explains in a satisfying manner the reason why, in every 
case known, the kinetic characteristics of the liquid 
(which can to first approximation be represented by its 
diffusivity) evolve with temperature in such a way as 
always to generate equilibration times of the order of 
experimental time scales before the thermodynamic 
crisis is upon us. 

Considerable headway on each of these problems was 
made by Gibbs and his coworkers [5,6,7] though they 
would be the first to recognize that the problems re- 
mained far from resolved by their work. Indeed, the 
"problem of glasses" remains a major challenge to the 
condensed matter theoretician and his experimental col- 
leagues at this time [8]. 

Before reviewing the contributions of Gibbs et al., let 
us consider briefly what is implied about the topology of 
the chemical potential hypersurface which must be rep- 
resentative of substances exhibiting the type of behavior 
described above. Here we merely rephrase much of 
what was written by Gibbs and his contemporary 
Goldstein in articles written now some 25 years ago 
[9, 10]. 

The fact that glasses are brittle solids at temperatures 
below their glass transition temperatures implies that the 
arrangement of particles taken up as a liquid cools can 
be described by a point in configuration space near the 
bottom of a potential energy minimum in this space 
[10]. If this were not so, the system would move in the 
direction dictated by the collective unbalanced force 
acting on it, and some sort of flow would occur. 
Notwithstanding the myth about medieval cathedral 
windows [11], this does not occur in glassy systems held 
at temperatures less than half their glass transition tem- 
peratures, even on geological time scales. On the other 
hand, the existence of the annealing phenomenon, in 
which the density and energy of a glass formed during 
steady cooling can change with time on holding at 
a temperature below but close to the "glass transition 
temperature" means that there is more than one such 
mechanically stable minimum available to the system. 
Indeed, there would appear to be an almost infinite 
number, of order e^, where N is the number of particles 
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in the system [12,13,14]. The minima obviously are 
distributed over a wide range of energies, the deeper 
minima usually being associated with higher densities. 
However, there are also many ways of organizing the 
same collection of particles into minima which differ 
negligibly in energy from one another. 

The fact that annealing proceeds more slowly the 
lower the temperature at which the annealing is carried 
out suggests that the process of finding deeper minima 
becomes more difficult statistically as the temperature 
is decreased. One arrives at the notion of an intercon- 
nected series of minima on a landscape of inconceivable 
complexity, in which increasing depth is coincident 
with decreasing population of minima, see Fig. 2 (from 
Ref. [2b]). The important implication of Kauzmann' s 
presentation in Fig. 1 is that for each system, at least for 
each "fragile" system, there must exist a statistically 
small number of minima at energies still well above that 
representing the crystal, and that these must set an abso- 
lute limit on the energy decrease achievable by anneal- 
ing the amorphous system. It is into one of these last 
few minima that the system tends to settle at the temper- 
ature where the excess entropy tends to vanish. The 
temperature characteristic of this ground state for amor- 
phous packing has become known as the Kauzmann 
temperature for that system, although the issue of the 
absolute value of the entropy, relative to that of the 
crystal that would be appropriate for a ground state 
system, remains unresolved. The issue is addressed in 
an important new paper by Speedy and Debenedetti 
[14] who succeed in evaluating quantitatively the num- 
ber and distribution of minima (each of which they term 
"a glass") for a model tetravalent system [15a, 15b]. 
This is an "inherent structure" [12] analysis which suc- 
ceeds in providing a quantitative description of the liq- 
uid thermodynamics in terms of the inherent structure, 
and which appears to identify the density of the "ideal 
glass" for the model considered. The authors argue that 
this ideal state would be reached by a second order 
transition during slow (but not necessarily infinitely 
slow [14]) densification, and conclude it would have an 
entropy in excess of the crystal at the same PV /NkT. 
Simple models of the two-state variety [16, 17] suggest 
the ideal glass would only be approached at K, despite 
linear extrapolations which would indicate otherwise 
[See Fig. 2(ii)]. The existence, in principle, of an ideal 
glass state has been disputed by Stillinger [12, 18]. 

The more minima per unit of energy, the larger the 
configurational component of the total heat capacity 
should be, hence the larger the drop in C p observed at 
the glass transition when ergodicity is broken during 
cooling. Of course, the change in C p reflects the density 
of minima at the level of the landscape at which the 
system gets trapped during cooling, and this level will 



depend on the height of energy barriers separating the 
minima as well as the total degeneracy, as will be 
described in the next section. Both features of the land- 
scape arise from the form of the interaction potential for 
the particles, but our knowledge of exactly how is in an 
elementary stage [12,13,14]. It is commonly found that 
fragile liquids, like the lactic acid of Fig. 1, have large 
changes in heat capacity at their T g , implying highly 
degenerate landscapes even quite close to 7^. 



(a) 
STRONG 



(b) 
FRAGILE 




Fig. 2. Sections through the 3N + I dimensional energy hypersur- 
faces of (a) strong and (b) fragile liquids. Z* is a collective configu- 
ration coordinate. Differences in the "density of minima" can be 
understood at an elementary level in terms of two-state models in 
which there are different increases in the number of distinct packings 
per elementary excitation event [16, 19] as represented in parts (ii), 
(iii), and (iv) of the figure for (ii) level of excitation, (iii) configura- 
tional heat capacity, and (iv) immediacy of the Kauzmann crisis, 
respectively, (from Ref. [2b] by permission). While this simple model 
clearly predicts, by extrapolation, an entropy crisis for fragile liquids 
at T > K it has been shown that a kinetic arrest about 20 % above 
T K does not satisfactorily resolve the Kauzmann paradox [16]. 



It is simpler to discuss the potential energy hypersurface 
as opposed to the chemical potential (free energy) hy- 
persurface, but it may also be less fruitful. For each 
interaction potential, there must exist a single im- 
mutable potential energy hypersurface for an TV particle 
system; however, the free energy hypersurface will de- 
pend on temperature. While a hard sphere system will 
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have a potential energy hypersurface that is totally de- 
generate — all configurations have zero energy — all con- 
figurations do not have the same free energy except at 
absolute zero. This is because the different configura- 
tions have different amounts of spare volume, and the 
free flight motion of spheres adjacent to such "loose 
spots" provides an entropy-generating mechanism. 
Thus different packings yield different pressures and 
different free energy minima. For the hard sphere sys- 
tem, these free energy minima have been evaluated by 
Dasgupta using a density functional approach [20]. 
Since we are most commonly interested in the behavior 
of systems at different temperatures, it seems that the 
most relevant hypersurface is the chemical potential hy- 
persurface. In any case, both of these hypersurfaces are 
quite impossible to conceptualize given that they exist in 
a space of dimensionality of the order of the number of 
particles. The attempt to represent them by two-dimen- 
sional slices, such as illustrated in Fig. 2, is a grotesque, 
but frequently practiced, oversimplification. It is permit- 
ted by the community only because it provides a way of 
thinking about such problems as the annealing of 
glasses and the configurational entropy of disordered 
systems. 

The individual minima on the free energy hypersur- 
face are the configurational microstates, or "configu- 
rons" [21] of the system. Note that the configurational 
entropy of the system is related to the number of minima 
accessible (in the thermodynamic sense) to it at a given 
temperature, irrespective of whether the system has time 
to "explore" them all. Thus a system that has become 
confined to a single minimum (which we call a "glass") 
still has an entropy in excess of the crystal, as shown in 
Fig. 1 by the bending over of the curves at the tempera- 
ture r g . This is problematical because statistically a sys- 
tem confined to a single state should have unit probabil- 
ity apart from vibrations. This was recently considered 
by Bowles and Speedy [22] who elaborate on the 
distinction to be made between statistical entropy and 
thermodynamic entropy, and conclude that it is the 
latter, indicated in Fig. 1 by the positive entropy values 
below r g , which is important in determining equilibria 
between physical states. 

While it is not obvious, the total number of config- 
urons per TV-particle system may be not very dependent 
on the nature of the interaction potential or even the 
organization of the particles into bound groups, i.e., 
molecules [15c]. Numbers roughly exponentiate with TV 
[12-14], e.g., exp(1.2A0 for the tetravalent model [14]. 
In this case the total entropy available by raising the 
temperature enough for all minima to be accessible, is 
£ B lnw, hence about £ B ln(expA0 i.e., about R entropy 
units per mole of particles. The "height" of the land- 
scape will then be different for liquids of different fra- 



gility, because of the different changes of heat capacity 
AC P they exhibit at T g . The quantity AC P is the heat 
capacity increment due to gaining access to the config- 
urational states, and it is its integral over the temperature 
range T K - T u which must amount to about R per mole of 
heavy atom centers. Thus T u , the temperature corre- 
sponding to the upper limit of the landscape (its 
"height"), will be high relative to T g for strong liquids, 
which have small values of AC P . Thus the liquid state 
"stretches out" with decreasing fragility, and crossovers 
to free diffusion behavior will occur at higher tempera- 
tures, as seen in recent studies [3, 4]. This matter is 
treated in more detail in coming publications. 

To be in equilibrium (except with respect to crystal- 
lization), a system must be able to visit, move between, 
a representative subset of the minima characterizing its 
chemical potential hypersurface. In this case, it is not 
correct to speak, as is often done, of a system well below 
its glass transition temperature as existing in any one of 
a manifold of "metastable states." These are 
"mechanically stable" but not "metastable" states. When 
held at a temperature near but below the glass transition 
temperature established during normal cooling, the sys- 
tem "anneals" by exploring the lower energy minima 
which were inaccessible time-wise during the initial 
cooling. Since, for most potentials, these are minima in 
which the particles are more densely packed, the glass 
volume will usually decrease during annealing. Even at 
constant density, however, annealing can occur since, in 
a complex system, configurations of different energy 
may have the same volume. Equilibration is made sim- 
pler by the fact that a macroscopic system consists of 
statistically independent nanoscopic regions such that 
the mean distance moved, by a single particle during 
equilibration, is only a fraction of a molecular diameter 
[23]. 

Annealing occurs more slowly at lower temperatures 
for two reasons. Firstly there are energy barriers to be 
crossed in the process of passing from minimum to 
minimum — the system collectively vibrates for increas- 
ingly long periods before some chance fluctuation 
(associated with exceptionally anharmonic excursions 
on the parts of some particles) permits it to rearrange, 
and hence to enter a new minimum. Secondly, there are 
entropic barriers because the lower energy minima are 
more distantly spaced, hence the number of successive 
rearrangements which must be made in order to arrive at 
one of the lower energy states must increase as temper- 
ature decreases. This aspect of relaxation has been 
considered recently by Mohanty et al. [24]. The combi- 
nation of barriers would seem to ensure that no matter 
how slowly a liquid is cooled (or equivalently how long 
a glass is annealed), the system will never reach the 
isoentropy condition in finite times (although T g itself 
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might appear to become independent of cooling rate 
[25]). Thus T K for the "configuron bath," like absolute 
zero for the "phonon bath," is inaccessible. 

This line of thought was made more quantitative by 
the development of the "entropy theory" of Adam and 
Gibbs [7], the usefulness of which we now discuss in 
some detail. 



3. Relaxation and Entropy 

At the time that the Adam-Gibbs equation was writ- 
ten down, it was very popular to explain the slow 
dynamics of supercooled liquids, and particularly of 
chain polymers, in terms of the free volume concept. 
The idea that a liquid's particles, or a polymer's seg- 
ments, move around at rates proportional to the amount 
of "elbow room" they find available, was both simple 
and satisfying, and the free volume concept was found to 
provide a convincing rationale of such important equa- 
tions as the WLF equation [26] which described the 
temperature dependence of relaxation times in a wide 
variety of polymeric liquids. The Adam-Gibbs theory, 
which was based on a modification of conventional tran- 
sition state theory to accommodate the notion that, in 
viscous liquids, the rearrangements over energy barriers 
must be cooperative and, further, that the size of the 
cooperating groups would necessarily increase with 



decreasing temperature, was running against the tide 
and was given little attention. Twenty-five years later, 
we see it in the ascendancy. 

The Adam-Gibbs approach led to an expression for 
the relaxation time which contains the excess (configu- 
rational) entropy, S c of the Kauzmann paradox, in the 
denominator of the exponent, a result of an inverse rela- 
tion between z*, the minimum size group of "beads" 
(rearrangeable units) capable of undergoing a rearrange- 
ment at temperature T, and the configurational entropy 
S c . The expression is 



Toexp hsrJ 



CD 



where A/jl is the conventional free energy barrier 
(per molecule in the cooperative group) to rearrange- 
ments, and C is a constant. The familiar departure from 
Arrhenius behavior comes from the temperature depen- 
dence of Sc which itself depends on the value of the 
configurational heat capacity. This is manifested at the 
glass transition by the change in heat capacity AC P . 

At the time of its publication, the Adam Gibbs 
approach was very appealing to the present author be- 
cause of the interesting behavior observed in some 
aqueous solutions under study at that time [27]. We 
encountered cases, such as that illustrated in Fig. 3, in 
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Fig. 3. Expansion coefficients for the liquid and glassy states of lithium acetate + 
water solution of mole ratio 1:10, as a function of temperature. Inset: a typical DSC 
trace obtained from a LiOAc ■ 10H2O sample showing the large jump in heat capacity 
at the same 7\ where the small decrease in a occurs. 
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which the glass transition temperature was associated 
with a zero or slightly negative change of expansion 
coefficient, Aa. Consistent with the second Davies- 
Jones relation [28], 



dT s /dp = VTAa/AC p , 



(2) 



these systems showed negative slopes of T g vs P, in 
striking contrast to the behavior of polymers and most 
simple liquids. Free volume concepts clearly do not 
make much sense for such systems, yet the same solu- 
tions showed large increases in heat capacity at T g (see 
Fig. 3), and hence behaved in a manner perfectly consis- 
tent with the Adam-Gibbs equation. Thus the Adam- 
Gibbs approach would appear to be more fundamental 
than the free volume approach since it applies irrespec- 
tive of whether or not entropy fluctuations are corre- 
lated directly or inversely (the cases of supercooled 
water, silicon, and Si0 2 , which all have density maxima 
below their T m ) with volume fluctuations [29]. 

We will not discuss the derivation of Eq. (1) except to 
note the general view of theoreticians that it cannot be 
understood (see, e.g., the contribution of E. A. Dimarzio 
in this volume). Notwithstanding this view, the equation 
appears to contain a lot of truth since it has excellent 
predictive capabilities. We will consider these capabili- 
ties under three conditions: a) relaxation very far from 
equilibrium, i.e., in the glass when the configurational 
energy is essentially constant; b) relaxation in the 
ergodic domain above T g ; and c) relaxation at tempera- 
tures near T g where the glass is annealing. 

3.1 Relaxation in the Non-Ergodic State 

For a fixed value of the configurational entropy, 
Eq. (1) predicts that such relaxations as may be ob- 
served, will have Arrhenius character, the slope of the 
Arrhenius plot being inversely proportional to the value 
of the entropy held constant. 

Evidence for the essential correctness of this predic- 
tion can be obtained from different sources. The first 
example is the Arrhenius variation of the electrical con- 
ductivity of ionic glasses. Well above T g , the ionic 
inverse conductivity follows the viscosity in its tendency 
to diverge near the Kauzmann temperature [30, 31]. 
However close to T g it tends to decouple, assuming a 
smaller temperature dependence, and then finally 
changes slope again at T g to assume its glassy state value 
[32]. As Eq. (1) predicts, the Arrhenius activation en- 
ergy is proportional to the amount of annealing which 
has been imposed on the glass since this lowers S c . The 
smaller S c becomes, the larger the glassy state activation 
energy. This correlation is more pronounced the more 
closely the conductivity is coupled to the viscosity [33]. 



A more direct example is offered by studies of vapor- 
quenched glasses of high excess entropy [34]. These can 
be studied at different entropy levels, fixed by annealing 
the deposits for different annealing times at higher tem- 
peratures. Equation (1) then predicts that, when the tem- 
perature is chosen to bring the relaxation into the exper- 
imental time window, the relaxation time (which can be 
obtained from sensitive calorimetric studies at the lower 
temperatures [28]) will be a linear function of the in- 
verse product TS C . The validity is demonstrated in Fig. 4 
[34]. As predicted, the slope is larger (and ris longer) 
the smaller the average value of S c . 

3.2 Relaxation in the Ergodic Domain 

3.2.1 Liquids Most of the tests of Eq. (1) have 
been carried out in the region of temperature above T s 
where the value of S c is an equilibrium quantity and one 
which changes systematically with the temperature. 
Eq. (1) predicts that under these circumstances, there 
will only be a single relaxation time vs temperature 
relation and that it will be a linear one for logr vs 
(TS c )~ l . An example is given in Fig. 5 for the case of 
tri alpha naphthyl benzene [35], for which the configu- 
rational entropy was obtained from differential scan- 
ning calorimetry. Log (viscosity) is seen to be linear in 
(TS C )~\ at low temperatures, and the curvature at high 
temperatures is as likely to be due to uncertainties in the 
assessment of S c [taken as either the difference between 
liquid and crystal entropies, ignoring differences in 
vibrational entropy at high T, (upper curve) or between 
liquid and glass at T g (lower curve)] as it is to failure of 
the theory. 

A more frequently used, though not as direct, test of 
the Adam-Gibbs equation is to develop the equation 
into a T Vogel-Fulcher-like form and then demonstrate 
that the T parameter has a value close to that of the 
Kauzmann temperature obtained by the purely thermo- 
dynamic route. In their original treatment, Adam and 
Gibbs made the simplest assumption for the excess heat 
capacity, which determines the configurational entropy 
temperature dependence, viz., that it is a constant. This 
then yields 



S c = AC p \nT/T K 



(3) 



Substitution into Eq. (1) then yielded the Vogel- 
Fulcher equation as an approximation, valid near T K . 
However AC P = constant does not describe many 
molecular systems. More accurate [36, 37] is AC P = 
KIT, from which 



AS(T) = K(T-T K )/TT K 



(4) 
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Fig. 4. Test of the Adam-Gibbs equation for the relaxation time of butyronitrile vapor- 
deposited samples (from Ref. [34] by permission). 
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Fig. 5. Test of the Adam-Gibbs equation for viscosity of tri-napthyl benzene at temperatures 
above T s . S c has been assessed in two different ways leading to two different plots (see text) , each 
of which is seen to be linear over a wide range of the variable (TS c )~ l (from McGill, Ref. [35] by 
permission). (A5 C in the paper of Magill is the S c of this paper.) 
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which leads to the Vogel-Fulcher equation as an identity 

T= T QXpDT /(T-To) = ToQXp(-F£) (5) 

where F is a fragility parameter, < F < 1 , and s - 
(T/T - 1), with T = T K . 

There has always been dispute concerning the validity 
of Eq. (5), and this has been revived with vigor recently 
in the light of a temperature derivative analysis by 
Stickel et al. [38, 39]. Stickel et al. show that, particu- 
larly for fragile liquids, there is a region at relatively high 
temperatures where the V-F equation fits quite well but 
yields a T that is considerably higher than T K and often 
is also larger than T g (a result which is unphysical and is 
associated with unphysical values of the r parameter). 
This analysis emphasizes the high temperature data 
whereas an analysis focusing on the last five to six 
decades in relaxation time (covering a small range of 
ordinate values on the Stickel plot) yields a lower value 
of To, one which usually agrees rather well with the 
Kauzmann temperature, while also yielding a physical 
(phonon-like) pre-exponential r . We will document this 
below for a large number of different glassformers. For 
intermediate liquids such as glycerol, even the Stickel 
analysis yields a T in good accord with T K [38]. 

In Table 1, we record T g , T K and T data for a com- 
pendium of liquids of different classes which have been 
studied over the years. The best fit values of T are 
reported for two different temperature regions relative 
to the glass transition temperature, where the data are 
available. When T depends on temperature, the value 
which should be compared with T K is the one associated 
with the most physical pre-exponent, e.g., r = 10" 14 s. 
This value is the inverse of the phonon frequency (/"(Hz) 
= 1/277T= 2 X 10 13 Hz) typifying the liquid quasi-lattice 
(which should set the attempt frequency for cooperative 
rearrangements). It is also the value of the pre-exponent 
found for relaxation in plastic crystals with very 
"strong" character (Arrhenius behavior) in which the 
extrapolation to \IT = is unambiguous [40]. 

The best fit r value, going with the value of T used 
in the final column comparison in Table 1 , is listed in the 
second final column so that the physical content of the 
T k /Tq comparison can be judged. If the number in this 
column is larger than 14, it suggests that the value of T 
used in the comparison with T K , although best-fitting 
the relaxation data, may be inappropriately low for the 
comparison. The values of the ratio T K /T would be 
unity if Eq. (1) properly described the liquid behaviour 
over wide temperature ranges, and if the basis for 
obtaining Eq. (5) from Eq. (2) is applicable. We note 
that in many cases of higher T s systems, AC P is not 
hyperbolic in T s , in which case Eq. (5) should only hold 
approximately near T s . 



Table 1 shows that T K /T values close to unity are 
obtained for glassformers with T g varying between 50 K 
and 1000 K. Since the liquids represented in Table 1 
range from molecular through covalent (Se, As 2 Se 3 ) to 
complex ionic oxides, we judge the case for the Adam- 
Gibbs approach to relaxation in glassforming liquids to 
be a strong one. 

An objection often raised to the Adam-Gibbs theory, 
and other approaches that similarly suggest that a ther- 
modynamic phase transition underlies the glass transi- 
tion, is that no experiment has detected a diverging 
length scale, £. However, Eq. (1), with S c ~l/z* ~l/(£*) 3 
together with Fig. 1, shows why a diverging length scale 
would be extremely difficult to detect. Since S c , even in 
the most fragile cases studied to date, rarely falls to less 
than one-third the entropy of fusion (Fig. 1, Refs. 1-30 
in Table 1), it must be expected that (z*)" 1 will only 
decrease to one-third of its initial value as T decreases 
from T m to T g . If £ (~z* 1/3 ) is of molecular dimensions, 
say 6 A at T m , then according to Adam and Gibbs, the 
characteristic length would only increase from 6 A to 6 
X (1/3 )" 3 , i.e., by ~2 A over this whole T m to T g range. 
Thus the Adam-Gibbs theory predicts that, even though 
£ should diverge as T — > T K , identifying a changing 
length scale by computer simulations, which can only 
probe much smaller decreases of S c than the above, will 
be an unrewarding endeavor — as indeed has been found 
[41]. 

3.2.2 Polymers Since many studies have demon- 
strated a relationship between the different canonical 
characteristics of relaxing liquids, and since aging 
(which is simply a slow and unwelcome approach to 
equilibrium at T < T s ) is a considerable problem in 
polymers, it is important to have some way of estimating 
how far above the ground state temperature a given 
polymer system is at ambient temperature. The problem 
with polymers is that reliable estimates of the Kauz- 
mann temperature are frequently not available because 
of the failure of the system to register a clean crystalliza- 
tion, hence a quantifiable entropy of fusion. However, if 
the Adam-Gibbs equation is as good as it appears from 
Table 1, then a reliable value of T would serve the 
purpose. We show how to achieve this in the following 
analysis. 

A great deal of relaxation data for polymers have been 
analyzed using the Williams-Landel-Ferry equation, 



log r/r g =Ci(T- T,)/[T - (T, - C 2 )] . 



(6) 



It is well-known that the Williams-Landel-Ferry and 
Vogel-Fulcher equations are mathematically equiva- 
lent. In the simplest case in which both are valid at all 
temperatures, values of T K would be available from the 
relation C 2 = T S - T assuming from before that T = T K . 



178 



Volume 102, Number 2, March- April 1997 

Journal of Research of the National Institute of Standards and Technology 



Table 1. Correlation of Kauzmann temperatures 


T K with Vogel-Fulcher T to values for various substances 








Substance 


r g a 


7k 


Ref. b 


T 
(high) 


T 
(low) 


Ref. 


T g /T Q 


T g /T K 


frag. 


(frag.)" 1 
D 


-logT d 

(-logTfa) 

C-logD ) 


T K /T Q 


1-butene 


58 a 


48 


1 


64(7?) 


[54(7?)] 


31 


1.07 


1.20 








>0.88 


2-methylpentane 


78 a 
83 


58 


2 


59 Cn) 
n-hex. 


60+5 Cn) 


32,33 


1.43 


1.38 


58 






0.97 


butyronitrile 


100 


81.2 


3 




58 (t d ) 


34 


1.72 


1.19 


47 


32 


16 


1.26-1 


ethanol 


95 


71 


4 b 




70-75 (t d ) 


35 


1.28 


1.33 






7 


1.0 




90 a 






80(t d ) 




36 


1.19 






2.7 


8.9 


0.93? 


methanol 


103 


64 ±5 


5,6 
9 


66(D) 


60±15(t d ) 


37(a) 
37(b) 


1.71 


1.61 




12.4 
12.4 


(7.1) 


1.06 
0.97 


n-propanol 


105 
100 a 


73 


7 


73.5(r,) 


73.5 (t d ) 
50.3 (t d ) 


38 
36 


1.42 


1.44 


33 
40 




11.7 
12.4 


1.00T 

1.45 


toluene 


126 


96 


8 


i03Cn) 


<108(t d ) 


39(a) 


>1.16 


1.19 




5.6 


(3.5) 


0.93 


tol + 17% benzCl 


126 






108(t d ) 


108 (t d ) 


39(b) 


>1.16 




107 




13.0 


-1.0 


ethylene glycol 


153 


115 


9 




109 (t d ) 


40 


1.43 


1.33 




16 


14.3 


1.05 


(ethan diol) 




119 


10 


125 c-n> 




41 


1.22 












1-3 prop, diol 


154 
145 a 


109 


10 










1.41 










1-2 prop diol 


172 


109 + 18? 


11 


i09(-n) 


114(t h ) 


42 


1.57 


1.41 


52 


17.8 


14.6 


1.11 


(T g =r g (l-3)+18 




(127) 






122(t d ) 


42 






52 


13.5 


13.2 


1.04 


glycerol 


193 
187 


135 


12 b 


121(t d ) 


128 (t h ) 
137 (t h ) 
127 (t d ) 
127 (t d ) 


42 

44 

38,9 

45 

45(b) 


1.51 
1.41 
1.52 


1.43 


53 
53 


19.5 

12.7 
8 


15.6 

14.6 
33 


1.04 

0.99 

1.07 

1.07-1 

1.11 


H 2 SO 4 *lH 2 


182 


142 


13 b 


146(a) 




46 


1.25 


1.28 








0.97 


H 2 S0 4 *2H 2 


169 


131 


13 b 


120(a) 




46 


1.41 


1.29 








1.09 


H 2 S0 4 -3H 2 


162 
155 a 


135 


13 b 


128(a) 




46 


1.27 


1.20 








1.05 


H 2 S0 4 *4H 2 


157 


133 


13 b 


136(-n) 
136(a) 




46 
46 


1.15 


1.18 








0.98 


triphen.phosfite 


205 


166 


14 


183 


186 


34 


1.12 


1.23 


160 


2.9 


13.3 


0.91 


PMS (disiloxane) 


165 a 


137 


15(a) 






47 




1.20 










Salol 


220 a 


167 


16 




135 (t h ) 


4d 


1.63 


1.31 




38.7 


24.5 


1.24-1 






157 


17,48 




141 (t d ) 


48 


1.56 


1.40 




33 


23.3 


1.1U 


orthoterphenyl 


244 


200 


18 




184 (t h ) 


49(a) 


1.33 


1.22 






17.7 


1.09-1 


and 










196 (ri) 


17, 32 


1.26 


1.26 


81 c 






1.04 


otp + 16% opp 










193 (r,) 


50 














dibutyl phthallatel 


179 




19 


151(-n) 


151 Cn) 


51 


1.18 








(3.6) 




(uncrystallizable) 










137 (t d ) 


32, 49(b) 


1.30 




69 




14 




m-toluidinc 


187 


151 


14, 8(b) 


153 (t d ) 


8(b) 






79 




13 


1.00 






157 


15(b) 




















propylene 


156 


125.8 


20 




130 (t d ) 


52 


1.20 


1.23 


104 


2.9 


13.1 


0.97 


carbonate 




128 


52b 




132.3 (t d ) 


36 












0.95 


Ca(N0 3 ) 2 *4H 2 


217 


200 
204 


21 

22 


205 Cn) 
201(o-) 


190(77) 


53 
54 

55 


1.14 
1.08 


1.09 








0.98 
1.05 
1.0 


Cd(N0 3 ) 2 *4H 2 


213 


198 


21 






56 




1.08 










fructose 


286 


210 


14 




206 (t e ) 


57 


1.39 


1.36 






13.5 


1.02 


glucose 


306 


271 a 


23 




231 
259(t d ) 


58 
43 


1.32 


1.13 
1.18 






12 


1.17 
1.05 
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Table 1. Correlation of Kauzmann temperatures 


r K with Vc 


igel-Fulcher T to values for various substances — Continued 






Substance 


1 $ 


T K 


Ref. b 


To 
(high) 


T 
(low) 


Ref. 


TJTo 


T,/T K 


frag. 


(frag.)" 1 
D 


-logTo 

(-logT/o) 

(-logDo) 


T K /To 


mannitol 


282 


236 


9,24 




















sorbitol 


266 


236? 


24 


212Cn) 




59 




1.13 


93 


8.6 




1.11 


(dulcitol) 




217 


24 




224 (t d ) 


9 


1.19 


1.23 




7.8 


14.3 


1.05 


sucrose 


323 


283 
287 


25 
14 




290 


58 


1.11 


1.14 
1.125 




0.154 




0.98 


trehalose 


388 




14 






57 










13.5 




phenolphtalein 


363 


310 


14 




274 (t e ) 


57 


1.32 


1.17 






13.5 


1.13 


selenium 


307 


240±10 


26 




25icn) 


32,63 


1.22 


1.28 


87 c 






1.04 


ZnCl 2 


380 


250 ±25 


27 


260(ti) 


180(tO 
236(t v ) 


60 

32,64 

32,65 


1.46 
1.61 


1.52 


30 c 
42.5 C 




14 
14 


0.96 
1.39 
1.06 


Li acetate 


401 


381 e 


28 


371(a) 




61 


1.08 


1.05 a 






14 


1.03 


As 2 S 3 


455 


265 


26 




237(T h ) 


44 


1.82 


1.93 






18.7 


l.OOl 


La 2 02B 2 3 


959 


845 


29 


864 


850 Cn) 


29 


1.12 


1.13 








0.99 


CaAl 2 Si 2 8 


1118 


815 


30 




805 (t|) 


62 


1.39 


1.37 








1.01 



a T g value based on the onset C p from adiabatic calorimetry, which is several degrees lower than scanning calorimetry or DTA-based values because 

of the much longer time scale. T g /T K values are based on 10 K/min DSC or DTA data for T g . 

b Indicates that the assessment of T K will not be found in the calorimetry paper cited, but rather in one of the authors articles, or students' thesis. 

(i) C. A. Angell and W. Sichina, Ann. N.Y. Acad. Sci. Vol. 279 (1976) p. 53. 

(ii) E. J. Sare, Ph. D. thesis, Purdue Univ. (1970). 

(iii) D. L. Smith, Ph. D. thesis, Purdue Univ. (1983). 
c Many values of m, defined as the slope of a Fig. 6 type plot at TJT - 1, are compiled in R. Bohmer, K. L. Ngai, C. A. Angell, and D. J. Plazek, 
J. Chem. Phys. 99 (5), 4201^209 (1993). Where the m value is used to obtain T Q via Ref. 32, the superscript c is attached to the m value. Such 
T Q values are associated with -logr of 14 by assignment. 

d -logTo is the value of To which is the best fit value for the T Q value cited. If -logr is numerically larger than the physical value of 14 (phonons), 
then the T Q value should be weighted up, and therfore T^/Tq should be weighted down. Where this is an important effect, the value of T^/Tq is tagged 
t or 1 to indicate the need for adjustment. For viscosity, the equivalent value of -Iog(i7 /P) is 3.5 and for diffusivity -log(D /m 2 s _1 ) is 7.55. For 
cases in which T Q is obtained from an m value via Ref. 32, the value of -logTo is 14 by assignment. 
e Unpublished data (Sichina and Angell) suggest this estimate of T K is too high, that AC P passes through a maximum and T K retreats to ~ 360 K. 

(T To value from conductivity. 

rj To value from viscosity measurements. 

t d T Q value from dielectric relaxation measurements. 

t e Tq value from tensile stress relaxation measurements and assignment t = 10~ 135 s. 

t h Tq value from ac heat capacity measurements. 

t l To value from longitudinal relaxation time from digital correlation spectroscopy. 

T h To value from Sherer-Hodge T $ analysis Ref. 44. 

t v To value from volume relaxation activation energy at T g and Ref. 32. 
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However, just as the Vogel-Fulcher equation usually 
does not provide an adequate fit of the data over the 
whole range of possible relaxation times, so does the 
WLF equation provide parameters which are frequently 
inapplicable outside the range of data taken, hence 
which are unreliable guides to the Kauzmann tempera- 
ture. As seen above (Table 1), for applications of the 
Vogel-Fulcher equation to molecular liquids, T is found 
to lie close to T K when the pre-exponent, r , of Eq. (1) 
lies close to 1CT 14 s, which is the physically meaningful 
value for barrier crossing (or de-trapping) processes. We 
have shown elsewhere [42] that for the case r = 1CT 14 s, 
C\ of the WLF equation should have the value 16 
because of the relationship (which is not generally to be 
found in textbooks), 



Ci = logTg/To. 



(7) 



For cases with r = 1CT 14 s, the value C\ = 16 follows 
from the definition r g = 10 2 s. Thus for cases in which 
the WLF parameters include Ci ~ 16, the ground state 
temperature T K can be obtained as T g - C 2 . Likewise the 
fragility can be taken as 



F = 1 - C 2 /T z = TJT Z = TJT Z 



(8) 



(where T is the Vogel temperature) since this provides 
a number which varies between zero and one with in- 
creasing fragility. The most fragile of the chain poly- 
mers, according to various studies [47] is polyvinyl chlo- 
ride (PVC) for which the long-established WLF 
parameters are Q = 16.2, C 2 = 25 K [48]. Since Ci has 
the physical value, T K can be obtained from C 2 and T g . 
The fragility accessed by the above expression, Eq. (7), 
is indeed close to unity, 1 - C 2 /Tg = 0.85. By contrast, 
polyisobutylene (PIB) has d = 16.5, C 2 = 104, and F = 
1 - C 2 ITg = 0.48. The general consistency of chain 
polymers with the strong/fragile pattern, and with to 
values of the order of the vibration period is shown by 
the selection of data in Fig. 6. 

3.3 Relaxation in the Non-Ergodic State Near T g 

Near T g , relaxation is complex because the quantity S c 
of Eq. (1) is time-dependent. Thus a measurement per- 
formed at a constant temperature will yield a relaxation 
time for recovery of the equilibrium state which is never 
a linear function of the displacement because S c itself 
relaxes according to an Adam-Gibbs equation. When S c 
finally reaches its equilibrium value, the relaxation time, 
of course, becomes time-independent. Behavior in this 
region has been treated in detail in the review by Hodge 
[49 ] and is additionally covered by the contribution of 
Hodge in this symposium issue. Here we note only that 



if the temperature of the isotherm falls below T K , then 
the relaxation time will necessarily diverge as waiting 
time increases. A diagram summarizing the relation 
between the accessible microstates, relaxation time and 
configurational entropy is given in Fig. 7. 
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Fig. 6. Tg-Scaled Arrhenius plots of segmental relaxation times for 
linear chain polymers based on mechanical relaxation, light scatter- 
ing, and NMR 13 C segmental relaxation time data. T g is defined by the 
temperature at which r = 10 2 s. The data extrapolations suggest an 
infinite temperature value of about 10" 14 s, consistent with Raman 
modes associated with short wavelength acoustic phonons which, in 
the absence of selection rules, dominate the Raman spectrum in this 
frequency range • [43]. The identifications of data sources on the 
diagram are as follows: INS = inelastic neutron scattering for the 
inverse wave vector, Q 1 = 1 A; [44] QUENS (Q Q ) = Quasielestic 
neutron scattering at peak of structure factor, Q ; [45] ILS = impulsive 
light scattering [46]. 



4. Systems With Discontinuous Changes 
in Configurational Energy and En- 
tropy 

In a minority of cases, mainly those with tetrahe- 
drally coordinated, inefficiently packed, ground config- 
urational states, discontinuous changes in configura- 
tional entropy, hence in relaxation time, seem to be 
possible. In the laboratory liquids so far discussed, this 
phenomenon appears to occur in the supercooled state 
in association with high crystallization rates, so the phe- 
nomenon is not easily studied. Furthermore, it is often 
not possible to be sure that the change is discontinuous 
(first order) rather than continuous through a smoothed 
higher order transition [50]. The clearest case which has 
been documented by computer simulation studies [51] 
as well as circumstantially by laboratory studies [52] is 
liquid silicon. The case of water is controversial 
[50,53,54,55,56,57]. Current studies [58] suggest that a 
similar phenomenon should occur in liquid beryllium 
fluoride at temperatures above the melting point (but 
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Fig. 7. Illustration of the relation between relaxation time, entropy, and excitation level on the potential energy hyper surf ace for a fragile 
glassforming system. Point 1 is in the free diffusion regime, unperturbed by any barriers to cross or traps to escape from. This is the regime of 
mode coupling theory validity. Around point 2, the melting point of this model glassformer, the system begins to "sense" the landscape, and 
it becomes increasingly enmeshed as T — » T $ at point 5. With relatively little excess entropy (the "lifeblood of the liquid state") remaining, the 
system falls out of equilibrium, becoming trapped in a single minimum (becomes a "glass") as its relaxation time rapidly increases beyond the 
normal experiment measurement timescales. The ideal glass is the configuration shich has energy within the lowest well of all (excluding crystal 
wells), which would be occupied at T K in a sufficiently slow cooling process. Then the configurational component of the total entropy would be 
* B lnl =0. 



also, unfortunately, above the normal vaporization tem- 
perature). It is hoped that, by control of pressure, some 
direct observations of the phase transition in the latter 
system may become possible. 

An important finding associated with these transi- 
tions, whether first order or continuous, is that the char- 
acter of the liquid changes abruptly on approach to the 
fully coordinated state and that the change is one from 
fragile to strong liquid behavior. This interesting phe- 
nomenology is in the early stages of its exploration. An 
account of the origin of the phenomenon in its different 
possible manifestations is given in the thermodynamic 
model of Poole et al. [59] and in the several microscopic 
models by Debenedetti and coworkers [60,61,62,63]. 



The entropy discontinuity, according to one parameter- 
ization of the Poole model, is shown in Fig. 8. 

There is an interesting possibility that the same phe- 
nomenology might be available in mesoscopic systems 
of the biopolymer variety. Evidence presented else- 
where [2c] suggests that when proteins unfold, there is 
a change in character from strong liquid to fragile 
liquid. Since in these cases there is no interference from 
crystallization phenomena, there is some hope that 
systematic studies of this provocative phenomenology 
might be made with such systems or with model 
systems based on them. Not only proteins behave in 
this manner but so do various RNA molecules. The 
attempt to correlate this aspect of biophysics (or, better, 
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Fig. 8. Discontinuity in entropy, and precursor effects at a liquid-liquid transition obtained at 
ambient pressure in one parameterization of the Poole model for water-like substances. This can 
occur when the hypersurface is characterized by the two megabasins [2c] separated by a 
substantial energy barrier. The megabasins usually seem to differ in topology, the high density 
one being characteristic of a fragile liquid, and the low density one, strong. (Adapted from Ref. 
[68].) Dotted line shows behavior for parameterization with stronger bonds. 



biopolymer physics) with the phenomenon of polyamor- 
phism in molecular systems, is in its infancy, but it bears 
directly on such important societal problems as "mad 
cow" disease, which involves refolding of "good" 
proteins into lethal polyamorphic forms [64,65,66]. This 
relation is discussed in more detail elsewhere [67]. 
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